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Abstract 

We develop the quantum inverse scattering method for the one-dimensional Hubbard 
model on the infinite interval at zero density. i?-matrix and monodromy matrix are obtained 
as limits from their known counterparts on the finite interval. The i?-matrix greatly simplifies 
in the considered limit. The new i?-matrix contains a submatrix which turns into the ra- 
tional i?-matrix of the XXX-chain by an appropriate reparametrization. The corresponding 
submatrix of the monodromy matrix thus provides a representation of the Y(su(2)) Yangian. 
From its quantum determinant we obtain an infinite series of mutually commuting Yangian 
invariant operators which includes the Hamiltonian. 
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1 Introduction 



Among the exactly solvable Id quantum systems the Hubbard model has probably the most 
interesting applications in solid state physics. Its elementary excitations, their dispersion 
relations and S-matrix at half filling |||] have been calculated exactly by use of the coordinate 
Bethe Ansatz in conjunction with the S0(4) symmetry of the model |2|, ^. The story of 
this development is long and originates in the seminal paper Q of Lieb and Wu. A recent 
overview is offered by the reprint volume 

Prom the point of view of coordinate Bethe Ansatz the one-dimensional Hubbard model 
appears similar to the fermionic nonlinear Schrodinger model. In fact, Lieb and Wu in 
their article used this analogy to obtain the Bethe Ansatz equations from Yang's earlier 
result ||6|, |^. Algebraically, however, the model seems to be more complicated. Nearly 20 
years passed before the basic tools of quantum inverse scattering method (QISM), i?-matrix 
and L-matrix were derived by Shastry |, g and by Olmedilla et al. |lO|, |ll|, |l2|, and it 
was shown only recently that the i?-matrix satisfies the Yang-Baxter equation (YBE) |13| j. 
The iZ-matrix and monodromy matrix of the Hubbard model have unusual features. The 
monodromy matrix is 4 x 4 rather than 3x3, as one might have guessed naively from the fact 
that there are two levels of Bethe Ansatz equations or from the analogy with the fermionic 
nonlinear Schrodinger model. It seems to be impossible to find a parametrization of the R- 
matrix, such that it becomes a function of the difference of the spectral parameters. For these 
reasons an algebraic Bethe Ansatz is difficult and was performed only recently by Ramos and 



Martins ||14|. 

There is another algebraic structure related to the Hubbard model on the infinite line. 



As was discovered by Uglov and Korepin [15| the Hubbard Hamiltonian commutes with two 
independent and mutually commuting representations of the Yangian Y(su(2)). Yangians are 



the quantum groups connected to rational solutions of the YBE 17, |18[. The Yangian 
invariance of the Hubbard Hamiltonian became likely after the observation that the S-matrix 
of elementary excitations at half filling is essentially a direct sum of two rational solutions of 
the YBE, each corresponding to the XXX spin chain [p!g|]. 

There is some hope that the Yangian symmetry might be used to obtain excitation spec- 
trum and n-point correlators of the Hubbard model in a way similar to the calculation of these 



quantities for the XXZ-chain by usage of its Ug(s/2) symmetry |20|. Such kind of approach 
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might also be applicable to an extended class of non-nearest-neighbour Hubbard models pl|], 



which have recently been shown to be Yangian symmetric, too [22|, and for which a QISM 
approach is unlikely to exist. 

The main concern of this Letter is to show how QISM and Yangian symmetry of the 
Hubbard model are connected. We benefit from the experience of one of the authors with 



the fermionic nonlinear Schrodinger model [p3| , 24 1. It turns out that the situation in case 
of the Hubbard model is to a large extent analogous. The Yangian symmetry reveals, when 
the model is considered on the infinite interval. Below i2-matrix and monodromy matrix 
are obtained as limits from their known counterparts on the finite interval. The i?-matrix 
greatly simplifies in the considered limit. The new i?-matrix contains a submatrix which 
turns into the rational i?-matrix of the XXX-chain by an appropriate reparametrization. 
The corresponding submatrix of the monodromy matrix thus provides a representation of the 
Y(su(2)) Yangian. This representation is identified as the Yangian representation constructed 
earlier by Uglov and Korepin using ad hoc methods. From the quantum determinant of the 
considered submatrix of the monodromy matrix we obtain an infinite series of mutually 
commuting Yangian invariant operators which is including the Hamiltonian. 

2 Quantum Inverse Scattering Method on the Infinite Interval 

The Hamiltonian of the one-dimensional Hubbard model is 



1\ / 1\ 1 



(2.1) 

— r \ ^ / \ ^ / 'i 

,1 3 

where cj^a and Cj ^ are annihilation and creation operators of electrons of spin a at site j 
of a Id lattice, and nj^j = c^j^Cja is the particle number operator. Since we want to study 
finite excitations over the zero density vacuum |0) of the infinite interval, we normalized the 
Hamiltonian such that H\{)) = 0. 

Starting point for the QISM for the Hubbard model is the exchange relation [|l^] 

n{\,fi)[Cm{\)(^sCM] = [£„(/i)0s A«(A)]7^(A,^), (2.2) 

where denotes the Grassmann direct product 

[A ®s B]^^^ps = (-l)[^(")+^(^)]^(^U,;3i?^5 (2.3) 
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with grading P(l) = ^^(4) = 0, P(2) = P(3) = 1. We adopt the expressions for the 



matrices TZ and £ in terms of two parametrizing functions a(A), 7(A) from ref. [11|. For later 
convenience, however, we shift the arguments of q;(A) and 7(A) by ^, such that we simply 
have a(A) = cos A, 7(A) = sin A. The L-matrix is 



•^m(A) 



(2.4) 



where /^^(A) = 7(A)(1 - n^ma) + ia(A)n 

ma ) 5»n(T 

(A) = a(A)(l - Umcj) - i7(A)nmCT, and h{X) 

is defined as 

smh2/i(A) _ [/ 

sin2A " T- ^^-^^ 

Due to space limitations we do not reproduce the i?-matrix here. It is 16 x 16 and contains 36 

nonvanishing entries, only ten of which are different modulo signs. The ten different entries 

are denoted by /Oj, z = 1, • • • , 10, in ref. [^]. They are rational functions of a(A), 7(A) and 

e'^(-'*). We provide a list and some basic formulae which have been used in our calculations in 



Appendix A. The (m — n)-site Hubbard model (m > n) is characterized by the monodromy 
matrix 

(A)/:„_2(A)---/:„(A). (2.6) 

It has been shown in ref. ^] that the logarithmic derivative of the graded trace of 7^n(A) at 
A = reproduces the Hamiltonian ( |2.1[ ) under periodic boundary conditions. Like Lm^X) the 
monodromy matrix 7^„,(A) satisfies ( |2.2D . In contrast to the classical case [25| a formulation 



of QISM on the infinite interval |26, ^] has so far only been possible for zero density of 
elementary particles. This is due to the complicated structure of the finite density (finite 
band filling) vacua formed by an infinite number of interacting particles. For the Hubbard 
model, there are four simple vacua, the empty band, the completely filled band and the half 
filled band with all spins up or all spins down. In the following we will consider the empty 
band |0) ( Cmo- 10) = 0) as reference state. Zero particle density means to consider only states 
with a finite number of particles in the empty band. Expectation values of the L-matrix with 
respect to this space have a finite limit for |m| 00, which formally can be obtained by 
setting normal ordered products of operators equal to zero ||2^, p7| , 

C{X) ^ V{X) = diag(-7(A)VW,a(A)7(A)e-'^(^),a(A)7(A)e-'*(^\-a(A)2e'^W). (2.7) 
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V{X) can be used to split off the asymptotics of 7^„(A). We expect the matrix 

%nn{X) = ViXy^'TmnWViXr (2. 

to have a finite hmit for m — > oo, n — > —oo. This hmit, 

r(A)= hm tnniX), (2. 



will be the monodromy matrix on the infinite interval. To derive an exchange relation for 
T(A), consider the asymptotics W{X,fi) oi Cm{X)®si^m{p) for large \m\ again by omitting all 
normal ordered products of Fermi operators. Then VF(A, ^)^'"(7^„(A) ®s %nn{fJ'))W {X, fj,)^ 
is expected to have a finite limit. 1^(A,^) is not just the tensor product V{X) (8s ^(/^)- Due 
to normal ordering there appear additional off-diagonal elements, 

W^(A,/x)i2,2i = W^(A,/u)i3,3i = -i7(A)7(/i), W^(A, ;u)i4,23 = -W^(A,/i)i4,32 = -i7(A)a(/i), 
W^(A, /x)24,42 = VF(A,^)34^43 = ia{X)a{^i), W{X, 11)23,41 = -W{X, 11)32,41 = -ia(A)7(/x), 
VF(A,/x)i4,4i = -e^W+'*('^). 

It follows from the exchange relation for the monodromy matrix that 

7^(A, fi)W{X, fi) = Win, A)7^(A, fi). (2.10) 

Using once more the exchange relation and the definition ( |2.8D of %nn{X) we obtain 

= [fmnil^) ®s %nn{X)]Un{n, X)'^n{X, /i)C/„(A, ^l), (2.11) 

where 

Um{X,fi) = iy(A,/z)-™[F(A)™ 0, y(/z)"^]. (2.12) 

Postponing the discussion of convergence for while we formally take the limits m, —n — > oo 
m (|ll|). Then by use of the definitions 

U±{X,^i)= liin ?7™(A,/i), ^W(A,;u) = C/±(;U,A)-l7^(A,;u)[/±(A,;u), (2.13) 

m— >±oo 

we arrive at the exchange relation for the monodromy matrix T(A) on the infinite interval, 
i?(+)(A,A^) [t(A) ®s t(/i)l = [t(/i) ®s t(A)l i?(-)(A,/i). (2.14) 
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The matrices U±{X,fi), their inverses and the i?-matrices iZ^^\X, p) can be calculated as 
functions of the pi's by utilizing the formulae provided in Appendix A| . The non-zero matrix 
elements of U±{X, fi) follow as 



U±{X, n)a/3,a/3 = f^±(A, ^)l4,41 
?7±(A,/i)i2,21 = f^±(A,/u)i3,3i = 
^±(A,^)23,41 = -C^±(A,//)32,41 



-P5 

P5- P4 



PW 



U±{X,fi) 



14,23 



-C/±(A,/U)i4,32 



1^6 



1/36 



P5 - Pa 

where pi = pi{X,p,). The corresponding matrices iZ^^^ are 
TZ^+\X,fi) = TZ^-HX,p) = 



U±{X,p)2i,A2 = C^±(A,^)34,43 



P3 - /5l 
^ iP2 
P9 ' 
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(2.15) 



The reader is urged to compare this expression with the i?-matrix on the finite interval |11]. 
Instead of the 36 nonvanishing elements of the original i?-matrix we have only 18 nonvanishing 
elements here, which brings about simpler commutation relations between elements of the 
monodromy matrix. 

From our experience with the finite interval case we know that the monodromy matrix is 
of the following block form, 



r(A) 



/ Z^ii(A) Cn(A) Ci2(A) Z^i2(A) \ 

Sii(A) An(A) yli2(A) Si2(A) 

B2l{X) A2l(A) ^22(A) B22[X) 

V ^21 (A) C2i(A) C722(A) D22{X) j 



(2.16) 
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where on the finite interval A(X) corresponds to the su(2) Lie algebra of rotations, D{X) 
corresponds to the ?7-pairing su(2) Lie algebra and the blocks B(X) and C(A) are connected 



to each other by particle- hole and gauge transformations [28|. Using the explicit form of the 
matrices iZ^^^ the exchange relation ( |2.16| ) implies 

P9(A /i)pio(A,M) lA^^x),A,si^^)]=A,^i^,)A^s{X)-A,piX)AM, (2.17) 

i.e. the commutation relations between the matrix elements of A{X) are decoupled from the 
rest of the algebra. This fact will be crucial for the derivation of the Y(su(2)) Yangian 
representation below. 

So far we avoided to comment on the analytic structure of the exchange relation ( p. 141) . 
This is indeed a delicate point. Along with the calculation of the limits Uj-{X,fi) we obtain 
convergence conditions. The convergence conditions for C/_|_(A, /i) and U-{X,fi) are comple- 
mentary. This situation is typical for the QISM on the infinite interval. Thus the first 
equation in ( 2.15| ) is only formal for the present time. We conjecture, however, that analytic 
continuation in A and /x respects the exchange relation ( |2.14[ ) with the possible exception of 
X = fi (modulo periods), where singular terms (like 5{X — fi) ) may destroy the first equality in 
( 2.15| ). We know from our experience with the fermionic nonlinear Schrodinger model |2^, 24] 
that such singular terms are irrelevant for the derivation of a Yangian representation from 
the exchange relation ( |2.14 ). 



3 Yangian Symmetry 



The Y(su(2)) Yangian |16, 17, 18] algebra is generated by six generators (n = 0,1; o 
1, 2, 3), satisfying the following relations; 



Qq: Qo 
Qq, Qi 

Qi,Q\ 



Qq, Qi 



+ 



Q1,Qi , Qq,Qi 



le' 



abc 



where k is a nonzero constant, a" (a = 1,2,3) are the Pauli matrices, /"'"^ 
antisymmetric tensor of structure constants of su(2), and ^"''^de/ _ jadk jbel jcfm jklm ^ Here 
and in the following we are using implicit summation over doubly occuring indices. The 



(3.1) 
(3.2) 

(3.3) 
is the 
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bracket { } denotes the symmetrized product 



i^^l J ■ ■ ■ ; Xmj — I / y ^crl ' ' ' X(jm ■ 

mn < ^ 



ml 



(3.4) 



cr£S„ 



Now we will show how to obtain a representation of Y(su(2)) from the submatrix ^(A) 
of the monodromy matrix ( 2.16| ). Introducing the reparametrization 



v{X) = -2i cot 2A cosh 2/i(A), 



the prefactor on the Ihs of (2.17) becomes 

P9(A,^)pio(A,Ai) 



P3(A,^)/34(A,m) - /'2(A,m)^ 
and we can write (2.17) in matrix form as 



v{X) - v{n) 



(3.5) 



(3.6) 



(if/ + {^(A) - v{,,)}V) [A{\) ® A{^l)] = [Aifi] ^(A)] (if/ + {v{X) - vif,)}V) . (3.7) 
Here "P is a 4 x 4 permutation matrix, (ill + {v{X) — v{fi)}'P) is the rational i?-matrix of 



the XXX spin chain. (3/7) implies that Aaf^{X) is a generating function of the Y(su(2)) 
Yangian. The Yangian generators Qq and Qf are the first coefficients in the asymptotic 
expansion [ 



2J 



oo -. / 3 N 

A^piX) = 1 + iC/ E ^Twm E Q>a(s + Q'Jap 

n=0 ^ ' \a=l J 



(3.8) 



In order to obtain compact expression for and Q\ we introduced the abbreviations = 
— (T^, CT^ = o"^, (7^ = (T'^.Q There are several possibilities to perform the limit v{X) — > oo. 
However, we found that only one of these yields finite results for Qq and Qf. We have 



to take Im(A) oo and further have to choose the proper branch of solution in eqn.(2.5) 
which determines h{X) as a function of A. Some of the details of the calculation are given in 



Appendix B . The final result is 



Qt 



(3.9) 



0</34,a4,7^*'7Ci,/3- (3-10) 



* The reader should not worry about this notation. It results from our choice of the L-matrix, which we 
took from ref. to facilitate comparison with earlier work. It is easy to introduce a slight change of the 
L-matrix, compatible with the exchange relation, such that a" is replaced by cr" in (B.8). 
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In our conventions the constant k in (3.3) is equal to ill. Comparing (3.9), (3.10) with the 



result of Uglov and Korepin |15], we find complete equivalence. 



Eo = Ql + iQl Fo = Ql- iQl, Ho = 2Ql (3.11) 
Ei = Q\ + iQl Fi=Q\- iQl Hi = 2QI (3.12) 

where the expressions on the Ihs are taken from the paper of Uglov and Korepin. 



Associated with the exchange relation ( ^.71) we can consider the quantum determinant |3C, 



DetgA(A) = Au{v{X))A22 {v{X) - iU) - Ai2{v{X))A2i {v{X) - iU) , (3.13) 

which is in the center of the Yangian and provides a generating function of mutually com- 
muting operators, 

[DetqA{X),Aap{fi)] = 0, [DetqA{X),BetqA{fi)] = 0. (3.14) 

The second of these equations is of course a consequence of the first one. Performing again 
the asymptotic expression in terms of v{X), 

oo ^ 

BetgAiX) = l + iUj2 ^TxY^^n, (3.15) 

n=0 ^ ' 

we obtain /q = 0, /i = iH, i.e. the Hamiltonian is among the commuting operators. All the 
conserved operators are Yangian invariant by construction. It will be interesting to investigate 



their relation to the formerly known conserved quantities |g, ^, 12, 32], which were obtained 
for the finite periodic model. 

In closing this section we shall add a comment. The Hubbard Hamiltonian on the infinite 
interval is invariant under the transformation 

The Yangian generators Qq and Qf, however, are transformed into a pair of generators Qq 
and Qf of a second, independent representation of Y(su(2)) ||T^. These two representations 
mutually commute. Therefore they can be combined to a direct sum Y(su(2))0Y(su(2)). 
The reason why we get only one of these representations from our QISM approach is that, in 
order to perform the passage to the infinite interval, we refer to the zero density vacuum |0). 
This vacuum has lower symmetry than the Hamiltonian. It is invariant under the su(2) Lie 
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algebra of rotations, but does not respect the //-pairing su(2) symmetry of the Hamiltonian. 
A fuhy su(2)©su(2) invariant vacuum would be the singlet ground state at half filling |1£]. 
It seems to be yet a formidable task to formulate the QISM with respect to this state. 



4 Concluding Remarks and Discussion 

We have developed the QISM for the Hubbard model on the infinite interval with respect 
to the zero density vacuum. The i?-matrix ( 2.15| ) thus obtained is greatly simplified in 



comparison with the i?-matrix of the finite periodic model. Particularly, it reveals a hidden 
rational structure, which arises from a certain combination ( p. 6] ) of the functions pi. This 
structure was discovered earlier by Ramos and Martins ||l^ as part of the exchange relation 
for the Hubbard model on the finite interval. Note, however, that our reparametrization 



which is essentially unique, differs from that given in |14|. A comparison is obstructed 
by the fact that the authors do not expose their parameters Oj . Along with the simplified R- 
matrix we obtained the aymptotic expansion ( |3.8[ ) of the submatrix A{X) of the monodromy 
matrix, which naturally provides a representation of Y(su(2)) and generates an infinite series 
of mutually commuting Yangian invariant operators including the Hamiltonian. 

There is a number of interesting open problems related to the QISM on the infinite 
interval. The analytic properties of the i?-matrices TZ^^'^ ( |2.15| ) and of the monodromy matrix 
( |2.9| ) deserve further investigations. Only the submatrix ^(A) of the monodromy matrix has 
a limit for Im(A) — > oo. All other matrix elements diverge. It is therefore not clear at the 
present stage of investigation how to obtain creation operators of elementary excitations that 
are compatible with the Yangian generators Qq and Q\. Creation operators are indispensable 
for the discussion of irreducible Yangian representations p3| , p4| . 



Another interesting task will be the construction of Dunkl operators [33| associated with 



the Yangian representation discussed in this Letter. Dunkl operators are building blocks of 



Yangian generators |33, 35, 23, 24]. They are useful for the investigation of eigenstates. For 



the one-dimensional Hubbard model, although several attempts [34, 35] have been made, no 
satisfactory Dunkl operator is known. 
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Appendix A Relations Between the Elements of the i?-Matrix 

In this appendix we collect functional relations among the elements of the i?-matrix, which 
have been used in the calculation of the matrix U±{X, fj,). We begin with the defining relations 
of the matrix elements, which are 

= e^«(A)a(A^)+e-'7(A)7(/^), (A.l) 



P4(A,/i) 
/32(A,/i) 
/39(A,/i) 
P2(A,/U) 

^io(A,/u) 



e'7(A)7(/u) + e~'a{X)a{^i), (A.2) 
= _e'«(A)^(^)+e-'7(A)a(/i), (A.3) 



: e'7(A)a(;u)-e-'a(A)7(;u), (A.4) 

/32(A,m) 

P3(A,^) _ e'a(A)Q!(^) - e"'7(A)7(^) 



P2(A,/i) a2(A)-72(/x) 

P5(A,/i) ^ -e'7(A)7(^) + e~'a(A)a(/x) 
P2(A,//) a2(A)-72(^) 
Pg{X, fi) e-''[e'a(A)7(A) - e-^a(/i)7(/i)] 



(A.5) 
(A.6) 
(A.7) 



P2(A,/i) a2(A)-72(/i) 
where h = h{X) + h{ii), I = h{X) — h{fj.). There are relations among the pi functions; 

Pi(A,)u)/)4(A,^) +/39(A,^)pio(A,/i) = p2{X,pf, (A. 8) 

P3(A,^)p5(A,/i) - P6(A,/i)^ = p2{X,pf, (A. 9) 

pi{X,p)p5{X,p) + p3{X,p)p4{X,p) = 2p2{X,pf. (A.IO) 



We found that there is a set of relations "dual" to ( |A.1| )-( |A?7 ). Introducing a transformation 
which keeps A unchanged and substitutes /^ + f for p, we get the following transformation 
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rules; 



P^ 

r J- 


— Pa Pa 




— /5i /3q jOin jOin 


pq 


P2 


P6 P2 




P6 P2 P6 P2 




P3 


P5 P5 




P6 _^ P2 




P2 


Pe' P2 


P6 


P2 P6 ' 




Explicitly, we obt 


ain the relations 










Pb{\p) - Pa{\p) 
P&{\p) 




—p^n!(\)'y(ii) -i- p~^'y( \)n'( lA 






P3{X, p) - pi{X, p) 
P6{>^,P) 




e^'y(X)a(n) - e~^a(Xh(u), 


(A.12) 




/5io(A,/i) 

P6(A,/U) 




-e^a{X)a(n) - e~''-f(Xh(u), 


(A.13) 




P9(A,/^) 
P6(A,/i) 




-eS(A)7(^)-e-'^a(A)a(/.), 


(A.14) 




P5{^,P) 




-e'^a(A)7(/i) - e~^-/{X)a{p) 


(A.15) 




P&{\p) 




a^X)-a^p) 




P2.{\P) 




-e'^7(A)Q(/x) - e-^a{X)-f{p) 


(A.16) 




Pe{X,p) 




a2(A)-a2(^) 




P2{\p) 




e-'[e'^a(A)7(A) +e-'^a(//)7(/i)] 


(A.17) 




P&{\p) 




a^{X)-a^{p) 



which are shown by direct calculation. The relations ( |A.8D -( |A.1C| ) are invariant under 

Appendix B Asymptotic Expansion of the Elements of the 
Monodromy Matrix 

We shall explain below details of expansion of the monodromy matrix T(A) in terms of 
v{X)~^ . Tmn{X) satisfies the recursion relation 

'Tm+l,n{X) = Crn{X)Tm,n{X), %n,m{X) = 1, (B.l) 

where 

CmiX) = ViX)-'^-'CmiX)ViXr 

( (icot A)"-T+"™i (icot A)"-Tc„|2^j^e^™PW 

-i(i cot A)"™T c^^^^Q-^rnp(X) ^^n^^ 

c"^ t(icot A)'"™i^e-i'"PW c"^ .c^i . , 

mt^ ' cosA mT ™-i sin A cos A 

i^t t 1^1- 2m\ I \\-n„|e2^ -imfc{A) 
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(B.2) 



-icm] (i cot A)"-i ^ei™p(^) -ic^^cmi ii^k cot''" A ^ 

™ I mj, sin A cos A '"IV / sin A 

(icot A)-"-T+"mj, (icot A)~"™Tc^^|^e^™fcW 

-i(icotA)-""TcJ^^^-^e-^"^'=W (i cot A)-"'"T-"™i j 

Here we introduced new functions 

gifc(A) ^ _g2h(A) gip{A) ^ _g-2h(A) (g_3) 

which we adopted from the recent analytic Bethe Ansatz for the Hubbard model by Yue and 
Deguchi |§. It follows from and (Qthat 

smfc(A) = --^ + — , smp(A) = --^ - — . (B.4) 

In the limit \m\ oo, the above matrix Cmi^) converges in the weak sense to the identity 
matrix. Solving (BTTJ) iteratively we obtain 'T{X) as 



t(A) = ■■■C^+i{X)C^{X)Cm-i{X)--- 

= l + ^(4(A)-l)+^(4(A)-l)(A(A)-l) + ---. (B.5) 

k k>l 

To expand Cm{X) in terms of ?j(A)~^, we consider the limit Im(A) oo and, to begin with, 
expand each function in terms of e'^'^. For e"^'*^'^^ there are two possible choices of branch, 



g-2MA) = _^ sin 2A ± y 1 + (^^ sin 2Xy . (B.6) 

To achieve convergence of the matrix elements T^/s (a, /3 = 2, 3) we have to take the lower 
sign in ( |B.(]| ). For this choice we get 

e-ipW = l{e2i^_2e^'^ + 0(e'5'^)}, JL = — {e^*^ + ©(e^i^)}. (B.7) 
U v{X) U 



Now the leading terms in the sums in ( |B.5| ) are of order e^'^, e*^'^, • • •. Thus, from the first 
two sums in ( |B.5| ), we get the expansion of the matrix ^(A) up to order e^'^. Then the last 



equation in ( B.7 ) yields the required expansion in {v{X)) up to second order 
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